An effective theory of Feshbach resonances and many-body properties of Fermi gases 
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For calculating low-energy properties of a dilute gas of atoms interacting via a Feshbach resonance, 
we develop an effective theory in which the parameters that enter are an atom-molecule coupling 
strength and the magnetic moment of the molecular resonance. We demonstrate that for resonances 
in the fermionic systems ^Li and ^''K that are under experimental investigation, the coupling is 
so strong that many-body effects are appreciable even when the resonance lies at an energy large 
compared with the Fermi energy. We calculate a number of many-body effects, including the effective 
mass and the lifetime of atomic quasiparticles in the gas. 
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Due to the enormous progress made in trapping and 
cooling atomic gases, it is now possible to create degen- 
erate Fermi gases 0, . The fact that interactions be- 
tween atoms can be varied essentially at will by exploiting 
Feshbach resonances makes these systems candidates for 
exploring experimentally strongly correlated Fermi sys- 
tems, in which the magnitude of the scattering length a 
is larger than the inter-particle separation. Many aspects 
of such systems have been considered, among them the 
equation of state '3| and the cross-over between the BCS 
state and a Bosc-Einstein condensate 0, . 

In this paper we begin by showing how, in the spirit of 
Landau Fermi- liquid theory and effective field theories, 
low energy Feshbach resonances may be characterized by 
a small number of phenomenological parameters. Re- 
lated ideas have been proposed in Refs. 0,0, and here 
we provide a compact quantitative formulation. To do 
this we give a field-theoretic account of a Feshbach reso- 
nance in a two-atom system which may readily be gener- 
alized to the many-body problem. Next we examine the 
conditions under which the Feshbach resonance induces 
strong correlations in a gas of fermionic atoms. Finally, 
we calculate a number of many-body properties. 

The first step in formulating an effective low-energy 
theory is to identify the relevant degrees of freedom. 
These we shall take to be two species {a =ti.) of atoms 
both of mass m, and a single low-lying molecular state 
whose energy relative to that of two stationary atoms in 
the open channel will be denoted by 2i'. It is important 
to remark that this molecular state consists partly of a 
"bare" molecule in a closed channel (or channels), to- 
gether with a dressing cloud of high-energy atom pairs in 
the open channel. Dressing of the molecular state by low- 
momentum pairs, which gives rise to threshold behavior, 
will be calculated explicitly. In addition, the matrix el- 
ement g describing coupling of atoms to the molecular 
state is different from the bare coupling. As an example, 
we calculate the width parameter AB which enters the 
standard expression a{B) = abg[l — AB/{B — Bq)] for 
the scattering length a as a function of magnetic field B 
near a Feshbach resonance. Here Obg is the background 
scattering length due to non-resonant processes. The 
resonant contribution to the effective atom-atom inter- 
action at zero energy is —g^ji^v). This in turn is equal 



to 47r?i^ares/w, where Orcs is the resonant contribution to 
the scattering length, from which it follows that 
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Here Afi — 2dv/dB is the difference in magnetic mo- 
ments between two atoms and the molecule. The result 
(P) is independent of the detailed microscopic model, pro- 
vided there is only one low-lying molecular state. 

We now demonstrate how the parameters 2v and g 
entering the low-energy theory are related to bare quan- 
tities for an explicit microscopic model, with the Hamil- 
tonian 
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Here E^'"'' = 2i^baro -l-p^/4m, where 2i/barc is the energy 
of a bare molecule with momentum zero measured with 
respect to the energy of a pair of atoms at rest in the 
open channel, efc = hF' /2m is the kinetic energy of an 
atom, .gbarc(fc) is the bare molecule-atom coupling ma- 
trix element, V is the non-resonant interaction between 
atoms, and V is the volume of the system. The quantity 
2i^bare cannot be measured directly because the energy 
of the molecule is affected by coupling to the open chan- 
nel, so we eliminate it in favor of the energy of the 
dressed molecule. The energy of a molecule at rest in a 
vacuum is given by the poles of the molecule propaga- 
tor Z?vac(^) = ^ 2i^barc ^ nvac(^)] "'^j whcrc in a com- 
pact matrix notation the vacuum molecule self energy has 

the form Hvac(w) = 5baic^vac(t^)ffvac('^ ) i. Fig. 3] and 

Gva.c{'^) = {lo — k'^/m)^^. The quantity (?vac(fc,'^) is the 
effective coupling between two atoms and a bare molecule 
in a vacuum, including effects of V. It is given explicitly 
by gva.c{k,z = k'^/m - iS) = (0m|//|(/)^), where l^m) is 
the bare molecular state and |(/)^) the scattering state 
for two atoms with relative momentum 2k and energy 
fc^/m in the center-of-mass frame interacting via V alone 
6]. We separate out the threshold behavior by writing 
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and 5vac,o = ffvac(0,0). The ^/— a; term is due to the 
threshold in the density of states for the open channel 
at w = 0. For positive energy it gives rise to damp- 
ing of the molecular resonance, while for w < it re- 
sults in the binding energy of the molecular state vary- 
ing as close to the resonance. We shall assume 
that frequencies of interest are small compared with the 
characteristic frequencies that enter the coupling ma- 
trix elements 5vac(fc,'^)- One then finds D^^^{lu) = 



uj/z 



-u;{4tt) ^, where 



nvac(O) - 

z~-^ — 1 — 9Anvac(0)/9w. Physically, z is the renor- 
malization factor for the molecular resonance when low- 
energy atom pairs in the open channel are neglected. 
This shows that the resonance energy, which is defined 
by Re[D-^l{2,^)] = 0, is given by 



2iy = Z[2zybarc + nvac(0)] 



(2) 



for 2z^baic + nvac(O) > 0. Thus, the energy of the reso- 
nance is shifted, and the difference between the magnetic 
moments of the two atoms in the open channel relative 
to that of the molecule is A/i — 2dv/dB = zA/ibarc, 
A/ibarc = '2dh'ha.ic/dB being the difference in magnetic 
moments between two atoms in the open channel and 
the bare molecular state. The energy of the molecular 
state is given for small B — Bq hy 21^ = A/i(i? — Bq). 
For negative cu one finds 1 — dUvs^d^) / duj = z^^ + 
g^jj^ gTO'^/2/(87r|a;|^/2). The singularity for lu 0^ re- 
flects the fact that the wave function of a very weakly 
bound molecular state consists mainly of low-energy pairs 
of atoms in the open channel, with only a very small 
closed-channel component 

According to standard theory, the resonant contri- 
bution Arcs to the scattering length is 47rajgj3/m = 

g2^,^0Z?vac(w = 0) = -5L,o/[2'^bare + H^ac (0)] 6]. With 

the help of Eq. this may be written as — zg^^^j. q/(2i/), 
in agreement with Eq. |^ if the effective coupling be- 
tween atoms and a dressed molecule is given by 
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The discussion so far applies equally well for bosons and 
fermions, apart from additional statistical factors if the 
two atomic species are identical. 

As an application, we consider in the rest of the paper 
a two-component Fermi gas. We assume that the popula- 
tions of the atomic species are equal and denote the total 
density by n. A dimensionless measure of the importance 
of the many-body effects due to the resonant interaction 
is fepjarcslj where fcp is the Fermi momentum, and we 
now estimate this quantity for systems of experimental 
interest. Feshbach resonances for fermionic atoms have 
been reported for ^'^K and ^Li. For ''"K, there is a Fes- 
hbach resonance for atoms in the states 1 9/2,— 9/2) and 



|9/2,-7/2) with parameters Bq ~ 201.5 G, AB « 8 G, 
and Obg — 164ao ^, The only other state in the 
ground-state manifold to which the open channel can 
couple by the central part of the interaction is 1 9/2, —9/2) 
and 1 7/2, —7/2). Diagonalization of the hyperfine Hamil- 
tonian yields A/^barc = 1.78/iB for B w 200 G whereas a 
coupled-channel calculation gives A/i = 0.118/iB for the 
molecule dressed only by high-energy pairs Q . From this, 
we obtain A/i/A/Zbaic — 0-067 for ^°K, which demon- 
strates that the renormalization of the molecule energy 
due to the coupling to the open channels is significant 
even away from threshold. With these parameters we find 

/cplarcsl — 15'T-]^2 ^f/i^ where ni2 is the density in units 
of 10^2 cm~'^ and ep = kp/2m. The dimensionless cou- 
pling is therefore quite strong for ^"^K. For ^Li we do not 
know the value of A/i , so as a first approximation we take 
A/ibarc ^ 2/iB IJ, [13 appropriate for the bare molecule. 
Using the parameters AB « 185 G and abg ~ 2160ao, we 

find /cpjarosl = 1-1 x ^0^n^2^^^F/tA:]are- At a frequency uj, 
the resonant contribution to the effective interaction be- 
tween atoms in the open channel is given by g'^ /{to — 2v). 
Since in calculations of many-body effects typical ener- 
gies are of order ep, this interaction may be treated as 
being independent of frequency provided 2^" 3> ep . Our 
discussion shows that even when this condition is satis- 
fied, systems can be strongly coupled with /cp|arcs| 3> 1. 

The next question we address is how the molecule 
properties are modified by the presence of other atoms. 
Such effects have recently been considered by Combescot 
|l2j |. The molecule propagator in a medium is given by 
= - S^-- + 2/i - n{p,a;) = Z)o(p,w)"' - 
n(p, lS) where n(p, w) is the many-body self energy 
and 2/i is the chemical potential of a molecule. As 
a simple approximation for 11 we take the same pro- 
cesses as in the calculation of Ilvac except that the 
atom propagator is that in the presence of a medium, 
rather than that in a vacuum [13. This yields 
n(p,t^) = 5LreGmb(P'^)5mb(?3,t^) where G|^^(p,k,w) = 

["F (Cp/2+k) + "•F (^p/2-k) - 1] / (Cp/2+k + ?p/2-k - - «<5) , 

with ^fe = ej; — /i, describes the propagation of a pair of 
atoms with total momentum p, energy w, and relative 
momentum 2k in the presence of the other atoms. Here 
np(a;) = [exp(/3a;) + 1]^-'^ with i3~^ = k^T. The many- 
body atom-molecule interaction in the ladder approxima- 
tion is given by g^h{p,uj) .gvac(w)[l - A.G^hTh^{u})]''^ 

with AGmb = a;) - Gvac(i^) and = + 2/i - 

/Am. Tbg is the vacuum T-matrix for the back- 
ground interaction V which for low energies is given 
by Tvac = 47rabg/m. One therefore has (7inb(p, = 
ffvac(w) + C'(fcFabg)- For clarity, we ignore such correction 
terms in the following discussion of the molecule proper- 
ties assuming fcp|abg| «C 1 and thus gmb(Piw) = (?vac('^). 
Expressing n(p, w) in terms of Ilvac we obtain 

zD-\p,uj) ^ LO-Ep-g''^Q{~[b)^\~n(p,uj). (4) 
Here Ep — p^ /Am + 2v — 2^ is the difference between the 
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molecule energy and the chemical potential when dress- 
ing by low-energy atoms is neglected. The term 



n(p,c^) 

5^ 



(2^ 



G[„'^(p,k,c.)~Gvac(fc,c:') (5) 



gives the contributions to the self energy coming from 
the presence of the atoms [l^ ■ In calculating fl we have 
replaced \/1.9vac(fc, (j-") by g = Vzgvac.o since contribu- 
tions from high momentum states, where the difference 
between g^i^dk^uj) and gvac.o is appreciable, are cut off 
by the Fermi functions. This approximation will be used 
in the rest of the paper. The importance of Eqs. Q 
and (O is that they show explicitly how molecular prop- 
erties at low energies may be expressed solely in terms 
of experimentally measurable quantities g and A/x, with- 
out arbitrary assumptions about the behavior of the bare 
coupling matrix elements for high momenta. 

We turn now to properties of atoms described by 
the propagator G{p,uj)~^ — Gq{p,u!)~^ — Il{p,uj) with 



Gq{p,uj) 



S,p. The atom self energy is 



given in the ladder approximation by icun 
/3"^Tr[r(K, q, q, zw„-|-iWm)Go(fc, iuJm)] where K = p-|-k, 
q = (p — k)/2 and the trace denotes a sum over 
ujm = (2m + l)7rfcBT and integration over k [T3| . 
The matrix T{p,Lu,k,k') obeys the equation r(p, w) = 

(2) 

Kff(p,w) + Vcs{p,uj)Gl^l{p,uj)T{p,uj) with T4ff(p,w) = 

V + gl^j.^Do{p,uj)gha.ic- After some manipulations S can 
be expressed as 



d^k k-p .o 
-^gmh(k + p, ^ ,uj) 



^ [?iF(e)Aa(fc, e)D{k + p, e + cj) 

-nB(e + t^)G(fc,e)An(k + p,e + cj)] (6) 

with nsix) — [exp(/3a;) — 1]^"'^ and Ai^{k,uj) — 
—2lmD{k,uj) and Aa{k,uj) = — 2ImG(fc,w) being the 
molecule and atom spectral functions respectively. Us- 
ing the low energy approximations Tvac = 47rabg/m and 
= gvac,o, we obtain 



5mb(p,k, w) 



1 



5vac,0 



d^q AGnib(P,q,^)rbg 

(2^)3 1- AG„,b(p,q,c^)Tbg 



(7) 



For T ^ 0, the self energy due to the background 
scattering V is given in the ladder approximation by 
Sbgfew) = [2fc|/(3TO7r)][fcFabg-he'(fcFabg)^] O. Since 
D is given by Eqs. and one sees from Eqs. © 
and O that the atom self energy may be expressed in 
terms of g and v (or equivalently A/i and Bq). 

From S one can calculate the properties of the atom- 
like excitations given by the poles of G. We begin by 
describing results for weak coupling, i.e. fcplaresl <SC 1- 
For 2v ^ ep and T = 0, one finds that the effective mass 
of an excitation at the Fermi surface is given by 



m 
m 



l + T^n+(211n2-3)^ 



■g'l^vfn^ 



40e2 



(8) 



and the wave function renormalization factor Z by 



Z-\p) = \ 



(9) 



The second terms in Eqs. ©-© come from the fre- 
quency dependence of the molecule propagator, while the 
final ones are identical with Galitskii's results for a dilute 
Fermi gas with short-range interactions parametrized by 
a scattering length a — abg-l-ares 1^ O'^^ neglects abg, 
the last two terms in the expression for m* have the same 
dependence on v, and the last term is ~ 2A:F|arcs|'^/eF 
times the second one. Thus for small g, the second 
term dominates over the last one while, for the reso- 
nances in ^Li and considered above, the last term 
is the more important at densities of experimental in- 
terest. For example, for the resonance, the ratio 
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-1/3 
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which demonstrates that the effective in- 



is 

teraction is equivalent to a contact interaction between 
atoms of strength g"^ /2v, provided that v ^ ep. Note 
that for the resonances in ^Li and '*°K considered in this 
paper, one can have A;F|a| > 1 even when v':^ e^- Thus, 
the gas can be strongly interacting due to the Feshbach 
resonance while the interaction still can be regarded as 
effectively instantaneous with a scattering length a, in 
agreement with the earlier order of magnitude estimates. 

The imaginary part of the atom self energy may be 
calculated in a similar manner. There is a contribution 
of order g^ coming from creation of real molecules. This 
is given for T = by Im[S](p, ^p)] — Sfeparest^^^^ep^^^ x 
{1 -p/\/128TOiy + [ep/(4i^) - l]{2mvf/^/p} for - 
\/2mep < p < VErniy + y/2meF and zero otherwise. For 
21^ ^ ep, this on-shell process is unimportant for atoms 
with p ^ 0{kp). For momenta outside the range for 
which the above process is allowed, damping is due to cre- 
ation of particle- hole pairs. One finds for \p~ kp\/ k-p <C 1 
and T ^ Tp to lowest non-trivial order in Obg and g that 



ImS(p,gp 
ep 



-(fcpa) 



(1 



kp ep 



(10) 



where again a — Obg -I- Ores. For T — 0, this result agrees 
with the expression derived by Galitskii for a gas of par- 
ticles interactin g vi a a short range interaction with scat- 
tering length a ^ . 

We turn now to stronger couplings, and we shall ne- 
glect flbg compared with arcs. The important ingredient 
in Eq. © is the molecule spectral function 



-2lmll{p,uj) 



[uj-Ep- ReU{p, cj)]2 + Imn(p, oj) 



(11) 



For calculating properties of the system at T < Tp, typ- 
ical energies of interest are of order ep. For such ener- 
gies, Imll is of order g'^mkp, and therefore for large v 
(/cp I Ores I ^ 1), the denominator in Eq. (|ll|l may be re- 
placed by (2i/)2, and one recovers the weak-coupling re- 
sult (|10|) . With decreasing molecule energy ly, the imag- 
inary term becomes increasingly important. In Fig. ^ 
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FIG. 1: — ImE(p, ^p)/eF as a function of feplarcsl for various 
momenta. The solid curves are numerical results based on 
Eqs. Q and lllll . the dashed curves are the results if the effect 
of the medium on the scattering amplitude is neglected, and 
the dotted curves in Fig.(b) are the Galitskii results [Ea. l|10|l ]. 
All these results become identical for feplarcsl ^ 1- 

we show numerical results for the quasiparticle damping 
rate -ImE(p, ^p) based on Eqs. ® and idJ, for T = 
and a coupling strength g appropriate for the *'^K reso- 
nance discussed above. For feplarcsl 3> 1, the damping 
rate saturates because the 11 terms in the denominator 
dominate. For comparison, we also show the results ob- 
tained by neglecting medium effects in the denomina- 
tor. In a Boltzmann equation approach, this amounts 
to assuming that the differential cross section for two- 



atom scattering with relative momentum k is given by 
its vacuum expression a'^^^/{l + fe^a^os), where effective 
range contributions have been neglected. As Fig. (a) 
shows, for large feplirosl the effect of the medium on the 
scattering amplitude increases the damping rate signif- 
icantly as compared to the Boltzmann results using the 
vacuum scattering rate. Remarkably, the scattering rates 
are substantially greater than one would predict if scat- 
tering cross sections were given by the unitarity limit for 
two particles in a vacuum. This is due primarily to the 
reduction by Pauli blocking of the magnitude of Imll , 
and hence also of the magnitude of the scattering ampli- 
tude. Figure Fig. (b) also shows that the results agree 
with those of Galitskii for fepl flics | ^ 1- 

The main results of this paper are, first, that we 
have shown how to formulate a theory of Feshbach reso- 
nances that involves only low-energy observables. One of 
these quantities is the magnetic moment of the dressed 
molecule, which is very different from that of the bare 
molecule in cases of experimental interest. We have 
also shown that for Fermi gases of experimental interest, 
many-body effects are strong even when the Feshbach 
resonance lies at energies well above the Fermi energy, 
and that for many purposes the interaction induced by 
the resonance may be regarded as instantaneous. As ap- 
plications of the theory, we calculated a number of many- 
body properties of a two-component Fermi gas. One im- 
portant result is that, for strong interactions, scatter- 
ing rates can exceed substantially predictions based on 
unitarity-limited two-body cross sections. 

We are very grateful to John Bohn for providing us 
with his result for A/i for ""^K. 
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